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Abstract. We compare the tau neutrino flux arising from the galaxy and the earth atmosphere for 10 3 < 
E/GeV < 10 11 . The intrinsic and oscillated tau neutrino fluxes from both sources are considered. We find 
that, for E > 10 3 GeV, the oscillated v T flux along the galactic plane dominates over the maximal intrinsic 
atmospheric v T flux, i.e., the flux along the horizontal direction. We also briefly comment on the prospects 
for observing these high-energy tau neutrinos. 

PACS. 95.85.Ry Neutrino, muon, pion, and other elementary particles; cosmic rays - 13.85.Tp Cosmic-ray 
interactions 



1 Introduction 

The Milky way is one of the nearby astrophysical sources 
producing high energy neutrinos, besides the familiar earth 
atmosphere [Q. Measurements of galactic neutrino and 
photon fluxes could provide information about the distri- 
bution of matter and cosmic rays in the galaxy. Further- 
more, the above flux is also a background for the search 
■ of more distant high energy neutrino sources such as the 
AGNs and the GRBs. In this talk, we shall focus on the 
flux of v r . It is clear that the observation of astrophysi- 
cal u T , with a flux comparable to the flux of v e and , 
directly confirms the neutrino oscillations. In order to ob- 
serve galactic v T flux, it is essential to study the atmo- 
spheric background. We shall focus on the energy range 
E v > 10 3 GeV. 

We first discuss the distinction between intrinsic and 
oscillated neutrino fluxes arising from the Milky way and 
the earth atmosphere. We then present results for galactic 
and atmospheric tau neutrino fluxes. Finally we comment 
on the prospects for observing galactic tau neutrinos. 



2 Intrinsic and oscillated neutrino fluxes 

It is well known that the relative flavor ratio for astrophys- 
ical neutrinos at the source is approximately cj)®^ : $J : 
<f>1 T =1:2:0. The commonly accepted astrophysical 
processes for producing electron and muon neutrinos are 
(7>p) + P ~ * 7r± + A where pion further decays into elec- 
tron and muon neutrinos with the ratio <\P V : =1:2. 



On the other hand, the production mechanism for the tau 
neutrino is (7,75) + p — > D s + X where D s further de- 
cays into tau neutrinos. The production cross section for 
D s meson is much smaller than that for -n^ for center-of- 
mass energy y/s ~ (1-10) GeV. Hence (j)° T is suppressed 
compared to 4>v c an d wt^ ■ 

Although v T flux is rather suppressed at the source, it 
is not negligible at the detector. As neutrinos propagate to 
the earth, the oscillation effect takes place and the relative 
neutrino flavor ratio changes as a result. Let us denote the 
neutrino flux reaching the earth as <\> Va . Then [2] 

& 

where P a p is a function of neutrino mixing matrix U a i 
which connects the neutrino mass eigenstate to the flavor 
eigenstate. The a, @ run over e, fi and r. Assuming a bi- 
maximal mixing for U a i, one obtains for vanishing 5 and 

013 

A/2 1/4 1/4 \ 
P a0 = 1/4 3/8 3/8 . (2) 
y 1/4 3/8 3/8 J 

Such a probability matrix implies (f> Ue : <j> v : (\> Vt = 1 : 
1:1. This flavor ratio is applicable to galactic neutrinos 
since these neutrinos propagate through a distance much 
greater than the neutrino oscillation length. On the other 
hand, this ratio is not applicable to the atmospheric neu- 
trinos for E v > 10 3 GeV considered here. At this energy, 
the oscillation length for the atmospheric neutrino is of the 
order 10 10 cm, which is much greater than even the earth 
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diameter. Hence the atmospheric v T flux for E v > 10 3 
GeV must be dominantly an intrinsic one. 



3 The galactic and atmospheric tau neutrino 
fluxes 

The total galactic tau neutrino flux consists of intrinsic 
and oscillated components. Both components follow from 
the collisions of primary cosmic-ray proton with interstel- 
lar medium proton. The density of interstellar medium 
proton is taken to be n p = 1 cm -3 along the galactic plane, 
while the primary cosmic ray spectrum, 4> p (E p ), is taken 
to be |U 



1.7 (Ep/GeV)- 27 
174 {E p /GeV)~ 3 



for E p < E , 
iovEp > E 0) 



(3) 



where Eq = 5 • 10 6 GeV and (j)p(Ep) is in units of cnr 2 
s -i gr -r Q e y-i, We assume directional isotropy in 4> p (E p ) 
for the above energy range. A more recent measurement 
of cosmic-ray flux spectrum between 2 • 10 5 GeV and 10 6 
GeV agrees with the (j> p (E p ) given by Eq. @ within a fac- 
tor of ~ 2 in this energy range [Sj . The intrinsic galactic 
tau neutrinos are produced by p + p — > D s + X, with the 
D s meson decays into a r lepton and a v T , while the r lep- 
ton further decays into the second v T with other particles. 
We calculate the D s production cross section by the fol- 
lowing two approaches: (i) the perturbative QCD (PQCD) 
and (ii) the quark-gluon string model (QGSM) [Hj. In the 
PQCD approach, we use the CTEQ5 parton distribution 
functions [7] and apply a K factor, K — 2, to account for 
the NLO corrections [Sj. With the D s production cross 
section determined, one can calculate the v T flux using 



AEp^p{Ep)f{E p 



1 



d<7 



j{Ep) 



pp-*v T +Y 

dE 



where f{E p ) = R/X P p(E p ) with \ pp (E p ) the pp interaction 
length and R a representative distance in the galaxy along 
the galactic plane. We take R to be ~ 10 kpc, where 1 pc 
~ 3T0 18 cm. We have focused on the intrinsic tau neutrino 
flux along the galactic plane just to obtain the maximal 
expected tau neutrino flux. The matter density decreases 
exponentially in the direction orthogonal to the galactic 
plane, therefore the amount of intrinsic tau neutrino flux 
decreases by approximately two orders of magnitude for 
the energy range of our interest. Another component of 
galactic tau neutrinos comes from the oscillation of galac- 
tic muon neutrinos. The flux of the latter can also be calcu- 
lated by Eq. |TJ} with v T replaced by Vp jSJ- In this case, n 
and K^- are the dominant intermediate states that decay 
into muon neutrinos. Since <j)® is few orders of magnitude 
greater than <fi®^ while (f>®^ is approximately one half of 
<f>® , we recover the flavor ratio </>° e : </>° : <\P Vt =1:2:0 
for the galactic neutrinos at the source. From the previ- 
ous section, we have <j) VT = (j)® /2. The flux tfi l/T can be 
parameterized as 



E^{E) = 



_ f 1.5 • 10" 5 (£/GeV)- 2 - 63 for E < E u 



Galactic-plane v x flux 
Atmospheric v t flux 
GZK v. flux 



1.1 cr' 



s 10 



10"' 
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(4) % 



9.5 • 10" 4 (.E/GcV)- 2 - 95 for E > E x , 



(5) 



Log,„(E/GeV) 

Fig. 1. Galactic-plane, horizontal atmospheric and GZK tau 
neutrino fluxes under the assumption of neutrino flavor oscil- 
lations. 



where E x = 4.7 ■ 10 5 GeV. The E<f>„ (E) is in units of 



Having discussed the calculation of galactic tau neu- 
trinos, we now turn to the atmospheric tau neutrinos. As 
said before, for E v > 10 3 GeV, atmospheric v T flux has 
only the intrinsic component. We have used the nonper- 
turbative QCD approach mentioned earlier to model the 
production of D s mesons in the pA interactions. We have 
used the (f>p(E p ) given by Eq. Q and the Z- moment de- 
scription for the calculation of intrinsic tau neutrino flux 
[TU| . We obtain the atmospheric tau neutrino flux by solv- 
ing a set of cascade equations |11U12| 

The results for galactic and atmospheric tau neutrino 
fluxes, along with the GZK 13 oscillated tau neutrino flux 
[T4] . are presented Fig. 1, where we have used the notation 
dV„ T /d(log 10 E) . From the figure, we note that the 
galactic plane oscillated v r flux dominates over the intrin- 
sic atmospheric v T flux for E < 5 • 10 7 GeV, whereas the 
GZK oscillated tau neutrino flux dominates for E > 5- 10 7 
GeV. Quantitatively, the atmospheric v T flux in the hor- 
izontal direction is ~ 5 times smaller than the galactic- 
plane v T flux. Furthermore, the downward atmospheric v r 
flux is factor of ~ 8 smaller than its horizontal counter- 
part. Let us recall here that particle physics aspects of 
intrinsic tau neutrino flux calculation presented here are 
empirically supported only up to y/s ~ TeV, which corre- 
sponds to E p - 10 6 GeV. 

Before we move on, it is instructive to compare the 
atmospheric v T and fluxes. We plot these two fluxes 
for E v > 10 3 GeV in Fig. 2. The atmospheric flux is 
taken from Ref. ^B]- The atmospheric Vp flux has two 
components: conventional and prompt components. The 
former component is due to decays of and K^, while 
the latter component is due to decays of charm hadrons. 
Although produced more copiously than charm hadrons, 
7r or K meson does not decay efficiently in the air at suf- 
ficiently high energy. This effect disfavors the resulting 
neutrino flux. Indeed, in Fig. 2, the two components of 
flux cross roughly at E v = 10 6 GeV rjHl- We also observe 
that the atmospheric Vp and v T fluxes are comparable for 
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The intrinsic atmospheric v &v flux 




Fig. 2. 

fluxes 



Log 10 (E/GeV) 

The comparison of downward atmospheric 



and v T 



In conclusion, we have presented our calculations of 
both the galactic and atmospheric tau neutrino fluxes for 
E v > 10 3 GeV. The former flux is shown to dominate over 
the latter. Such a dominance is unique to tau neutrinos. 
The event rate for galactic tau neutrinos by observing the 
double showers (with 10 6 < E„/GeV < 10 7 ) is rather sup- 
pressed. Therefore, to observe the galactic tau neutrinos, 
it is desirable to develop techniques for identifying tau 
neutrinos at lower energies. 

H.A. and K.C. are supported in part by the Physics 
Division of National Center for Theoretical Sciences un- 
der a grant from the National Science Council of Taiwan. 
G.L.L. and J.J.T. are supported by the National Science 
Council of Taiwan under the grant numbers NSC91-2112- 
M-009-019 and NSC91-2112-M-001-024. 



E v > 10 6 GeV, since both fluxes are due to charm pro- 
duction in this energy range. 



4 The prospects of observations 

For downward going or near horizontal high-energy tau 
neutrinos, the deep inelastic neutrino nucleon scattering, 
occurring near or inside the detector, produces two show- 
ers |16| . The first shower is due to a charged current neu- 
trino nucleon deep inelastic scattering, whereas the second 
shower is due to the (hadronic) decay of the associated 
r lepton produced in the first shower. It might be pos- 
sible for the proposed large neutrino telescopes such as 
the IceCube to constrain the two showers simultaneously 
for 10 6 < E/GeV < 10 7 , depending on the achievable 
shower separation capabilities f^Jj. Here, the two showers 
develop mainly in ice. With such a detection strategy, we 
estimate the event rate for observing the galactic tau neu- 
trinos. In the above energy range for the tau neutrinos, the 
event rate in 1 km 3 water/ice Cherenkov detector is rather 
small, about ~ 5-10 -3 yr _1 sr _1 . Such a low event rate im- 
plies that one can only search for galactic tau neutrinos in 
the lower energy. Alternatively, one should consider other 
detection strategies. Since Fig. ^ shows that galactic tau 
neutrino flux dominates over its atmospheric counterpart 
for E v > 10 3 GeV, it is desirable to develop strategies for 
identifying tau neutrinos at TeV energies or even lower 
because of relatively large absolute v T flux. The present 
AMANDA search for all flavor neutrino-induced cascades 
is not tight enough to constrain/observe our predicted <j> v 

m 

We point out that the persistent dominance of galactic 
tau neutrino flux over its atmospheric background is a 
unique phenomenon among all neutrino flavors. Such a 
dominance does not occur for v e and v^. For example, the 
result in Sec. 2 tells us that (f)^ = <j> Vr for the galactic 
neutrinos. On the other hand, the atmospheric flux is 
much greater than the atmospheric v T flux for E u < 10 5 
GeV, as can be seen from Fig. 2. Hence the galactic 
flux no longer dominates over the atmospheric flux for 
E v < 10 5 GeV. 
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